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In the last years there have been measurements of energy independent modulations in underdoped 
cuprates via scanning tunelling microscopy. These modulations of around 4 lattice constants occur 
near vortex cores in samples with and without superconductivity, and in the pseudogap regime 
of some superconducting samples. Recently, it was proposed that the origin of the modulations is 
electronic charge ordered phases intrinsic to doped Mott insulators. We show that even if the system 
has not charge ordered, proximity to such an order will induce modulations at the same ordering 
wavevector. Such modulations are a signature of the real or incipient order. We construct ground 
state wavefunctions that include strong fluctuations of the competing order without having ordered. 
They exhibit deep roton minima at wavevectors of the reciprocal lattice of the incipient order due 
to proximity to the phase transition into the charge ordered state. Such minima, in turn, produce 
strong modulations. Moreover, the effects are generic and are present whenever an electronic system 
is near a crystalline phase. This leads to experimental signatures even if the electronic system has 
not solidified. The signatures will be ubiquitous near the dopings that stabilize the electronic order. 
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In the present note we would like to point out that 
the recent proposed explanation^ for STM modula- 
tions at certain doping dependent wavevectors in samples 
with and without superconductivity, and around vortex 
coresQ is of a quite general nature. The modulations 
around the vortex core were predicted in the context of 
the SO(5) picture of high T c |3], but their nature is more 
general than the specific model. We show more directly 
how the physics comes about by means of a quantum me- 
chanical density functional formalism. The proposal^ 
is of an enhanced charge susceptibility due to charge 
order or to a deep roton minimum at certain wavevec- 
tors. We show that roton induced modulations are always 
present near a solidification quantum transition regard- 
less of whether the electron system has charge ordered 
or not. Therefore the experimental signatures should be 
quite generic and rather common as long as one is close 
to an electronic crystal phase. 

That roton physics follows on quite general grounds 
can be inferred from the experimental properties of 
bosonic systems. Helium is a superfluid at zero temper- 
ature, but it is an almost solid barely melted by quan- 
tum fluctuations. An indication of its proximity to the 
solid phase is the deep roton minimum in the excitation 
spectrum0,0 or, equivalently, the peak in the structure 
factor at the same wavevector. On the other hand, ar- 
tificially engineered superfluids such as atomic BECsjfJ, 
which are not near a solid transition by experimental de- 
sign, do not have a peak in the structure factorQ and 
hence they will not have a roton excitation unless they 
are tuned near a Mott transition. 

In the case of the cuprates there is a plethora of 
possible Coulombic stabilized insulating phases: the 
low doping antiferromagnetic insulator, stripes 8], hole 
Wigner crystals and even a semimetallic orbital 



antifcrromagnet |9J , among others. Whenever one of these 
phases with charge order (stripes, Wigner crystals, charge 
density waves) becomes competitive in energy with the 
superconductivity, an incipient Bragg peak will start 
forming, producing enhanced fluctuations at appropri- 
ate wavevectors and thus a spike in the structure factor. 
This spike leads to a roton minimum in the excitation 
spectrum 0, In fact, the magnetic neutron (%, ir) 

mode(llT| is such a roton minimum in the spin rather than 
charge density excitation spectrum. It is due to proxim- 
ity to an electron crystal phase, the antiferromagnetic 
insulator, at very low dopings 9] . 
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describes quantum many particle systems [4j. v ■ pv/2 
is the kinetic energy operator, U[p] is the potential en- 
ergy operator, which in general can be a functional of the 
density operator. The density operator in first quantized 
notation is 



(2) 



with ri being the position of the ith particle. The velocity 
operator is in first quantized notation 
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with commutation relation 
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The ground state of the Hamiltonian Q could be a quan- 
tum liquid or solid depending on the interaction. Since 
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our main concern is how a crystallization transition is 
approached, we will suppose it to be a liquid, of course 
a superconductor for the case of the cuprates, but our 
considerations are more general. 

We study the nature of the density excitation spectrum 
of the fluid on quite general grounds by the methods 
of Landau and the Russian school (l2T|. We suppose the 
density to have a well defined average which we take to 
be uniform for simplicity. This will not change the nature 
of the considerations, although in real life the density can 
be modulated by the lattice. The average velocity is zero 
as we consider a system at rest. We Fourier expand the 
density and velocity operators about their averages: 
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where we are using the lattice normalization, the number 
of sites N instead of the volume V. Quantum mechani- 
cally, the velocity is proportional to the gradient of the 
local phase variable of the particle, v(f) = h/mVO, or in 
Fourier components 
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The density velocity commutation relation implies the 
well known density phase commutation relation 
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i.e. is the momentum conjugate to the density fluc- 

tuation pr. The Hamiltonian Q thus becomes 
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and spring constant by 
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The density energy excitation spectrum of the quantum 
liquid is 



\P%\ 
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where U? is the Fourier transform of the interaction and 
(fi- 2 fc 2 )/(4mpo)|Pfc| 2 is the elastic energy associated with 
changing the density. This last term comes from acting 
the Hamiltonian on the wavefunction and making sure 
that one keeps careful track of both density and phase 
degrees of freedom 13] in the full quantum mechanical 
kinetic energy term. 

We see that the quantum liquid is a collection of har- 
monic oscillators in momentum space for the density fluc- 
tuations. The mass of the oscillators is given by 
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The ground state energy is given by U(po) + ^£ huj^/2. 
We note that the Virial theorem implies that 
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where the structure factor is defined by 
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which is, of course, the Fourier transform of the density- 
density correlation function. For the quantum liquid we 
thus have 
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The ground state energy will have the density oscillators 
unexcited. Hence, if the ground state wavefunction of the 
electronic system without density correlations is \ipGs)y 
the ground state wavefunction with density correlations 
is 



= exp 



poh 2 k 2 
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The factor in front appears because it is the ground state 
wavefunction of harmonic oscillators. The ground state 
wavefunction has the oscillators shaking as little as pos- 
sible. We see that the factor suppresses density fluctua- 
tions. The wavefunction is essentially exact to quadratic 
order in the density fluctuations. This way of includ- 
ing correlations was pioneered by Bohm and PinesfTi) 
and is known colloquially as the RPA. This correlated 
ground state wavefunction is a Jastrow type wavefunc- 
tion analogous to the Laughlin Fractional Quantum Hall 
wavefunction^^ . 

One could play the same game as for the ground state 
wavefunction 117(1 and write down the excited state wave- 
functions for the density oscillators. On the other hand, 
these collective density oscillations are mostly not good 
approximations to eigenstates of the system since they 
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will readily decay into particle-hole pairs except at the 
longest wavelengths. Such long wavelength oscillations 
are plasmons, for which the present consideration leads 
readily to cj§ = Aire 2 po/m, as for the electron Coulomb 
interactions we have £/g = 4-7re 2 /fc 2 which dominate as 
k -> 0. 

The ground state wavefunction i|17|) does not contain 
any electronic crystalline order as (pg) = for all k. 
When the electron system solidifies, either spontaneously 
or lattice induced, it will exhibit Bragg peaks at the recip- 
rocal lattice vectors corresponding to the electron crys- 
tal. These peaks mean that (p^) ^ 0. We are interested 
in the case when the electron system solidifies sponta- 
neously, and in particular, in its behavior as it approaches 
such a point. 



of the interaction we get the structure factor and the 
density excitation spectrum, which we plot in Figures 1 
and 2, where the roton wavevector was chosen arbitrarily 
as 1, a was chosen as 1/2, the roton gap A roton as 1/4 
and all other dimensionful constants were chosen to have 
unit scale. 

We consider the case when an external potential — Vjj 
acts on the electrons thus adding a term 



Tie 



k 



(19) 



The external potential might be the ionic lattice or it 
might be the impurity potential. The effect of the extra 
term in the ground state wavefunction is to move the 
equilibrium position of the density oscillators from (p-A = 
to 
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FIG. 1: Structure factor vs. momentum. The roton maxi- 
mum is easily visible. 

As the electron fluid approaches spontaneous crystal- 
lization with reciprocal wavevector Q for the incipient 
order, there will be enhanced density fluctuations dp 2 ?!) 
at such a wavevector. This follows because proximity to 
the crystallization softens the density oscillators of the 
electron liquid, —> 0, which makes uj^ — » 0. Hence 

the interaction Ur has a minimum at Q. The softening 
of the density mode thus produces a peak in the struc- 
ture factor, S(k), that diverges as one approaches the 
point when the system solidifies. This peak will drag the 
spectrum down according to the equation l(T4*|l creating a 
roton minimum [jl Il0j| that can be thought of as an exci- 
ton. We thus have as a simple phenomenological model 
that captures the physics 
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where nonzero a makes a roton minimum, and A ro ton is 
the roton gap, which will collapse at the crystallization 
transition as the smallness of the gap of such a roton is 
a measure of proximity to crystallization. From the form 
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We point out that when V? is the lattice potential, this 
is just the electron fluid acquiring the periodicity of the 
underlying lattice with appropriate Bragg peaks. 




FIG. 2: Density oscillation frequency vs. momentum. The 
plasmon gap and the roton minimum are easily visible. 

When the electron fluid is near a solidification tran- 
sition with reciprocal lattice vector Q for the incipient 
order, there will be a large enhancement of the density 
component at or near wavevector Q due to the roton 
minimum unless Vt = for the region around the min- 
imum. In figure 3 we plot the enhancement of the den- 
sity Fourier component for a constant potential which 
could correspond to a localized impurity. The scale of 
the potential in the figure was chosen arbitrarily as 1. 
On the other hand, we stress that the effect is general 
and the enhancement will be there for the lattice po- 
tential. Quasiparticles will scatter strongly of the roton 
enhanced density Fourier components at and around Q. 
These corroborates the previous suggestions hi and we 
thus identify this enhancement with the STM|2j bias in- 
dependent signatures. 
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We point out that our result shows the generality of 
the roton induced modulations near an electronic crys- 
tal phase effect, and that charge order is not necessary 
for the observation of the experimental signatures. The 
physics follows because it is very easy to make density 
fluctuations at the roton wavevectors. Since the elec- 
tron density is very high at roton wavevectors, there will 
be an enhanced charge susceptibility at such wavevec- 
tors leading to modulations [lj as long as the system is 
near enough to an electronic crystal phase. The phe- 
nomenon can be observed in superconducting samples 
which might not have the charge order. For example, in 
underdoped samples if superconductivity is lost at a tem- 
perature not far above the roton gap, these modulations 
will be present. Since T c collapses with underdoping, the 
modulations will be easier to observe for small doping. 
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FIG. 3: Density Fourier component vs. wavevector. The en- 
hancements of the density Fourier components at and around 
the roton wavevector are easily visible. 

Measurements of metallic like conductivity below 
the Neel temperature in underdoped cuprates, Hall 
effect [TU transport measurements, and the different 

photoemission properties of nodal vs. antinodal 

quasiparticles[16j have as the simplest explanation that 
there are two types of electronic carriers in underdoped 
cuprates. It is certainly suggestive that the experimen- 
tal wavevectors for the modulations are in the antinodal 
directions 0- The roton model will be a two fluid system 
as there will be low energy quasiparticles near the roton 
minimum which will behave more like excitations of the 
competing insulating order, while low energy quasiparti- 
cles coming from other wavevectors will tend to behave 
like superconducting Bogolyubov-BCS excitations. 

We finally point out that there are reports of a 
Boomerang effect in the cuprates 17]. This consists in T c 
becoming proportional to superfluid density, n s , in the 
overdoped regime as a consequence of a decreasing num- 
ber of superconducting electrons. This suggests a loss of 
spectral weight of the superconducting phase to a com- 
peting order. If the roton minima are indeed signatures 
of incipient electronic crystal phases and the competing 
phase on the overdoped side is an electronic crystal, there 



will be similar modulations in the overdoped regime ob- 
servable in STM experiments. 

We would like to mention that recently there was also 
a proposal that the STM modulations are a Copper pair 
crystal[l8|. On similar physical grounds as expressed 
in|l|, we believe such a crystal would not occur, but we 
emphasize that the two proposals are differentiable ex- 
perimentally according to at which doping one observes 
specific modulations. We also stress that the physics we 
have described is more general than either proposal be- 
cause as long as there is an incipient solid order, there 
will be roton minima at its reciprocal wavevectors and 
the concomitant modulations. 
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